Quantum Walk on a Line with Two Entangled Particles 
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We introduce the concept of a quantum walk with two particles and study it for the case of a 
discrete time walk on a line. A quantum walk with more than one particle may contain entangle- 
ment, thus offering a resource unavailable in the classical scenario and which can present interesting 
advantages. In this work, we show how the entanglement and the relative phase between the states 
describing the coin degree of freedom of each particle will influence the evolution of the quantum 
walk. In particular, the probability to find at least one particle in a certain position after N steps 
of the walk, as well as the average distance between the two particles, can be larger or smaller than 
the case of two unentangled particles, depending on the initial conditions we choose. This resource 
can then be tuned according to our needs, in particular to enhance a given application (algorithmic 
or other) based on a quantum walk. Experimental implementations are briefly discussed. 

PACS numbers: 03.67.-a, 03.65.-w, 05.30.-d 



Quantum walks, the quantum version of random walks, 
were first introduced in 1993 Q and have since then been 
a topic of research within the context of quantum infor- 
mation and computation (for an introduction, see 2]). 
Given the superposition principle of Quantum Mechan- 
ics, quantum walks allow for coherent superpositions of 
classical random walks and, due to interference effects, 
can exhibit different features and offer advantages when 
compared to the classical case. In particular, for a quan- 
tum walk on a line, the variance after N steps is propor- 
tional to N, rather than y/~N as in the classical case (see 
Fig. ^) . Recently, several quantum algorithms with op- 
timal efficiency were proposed based on quantum walks 
, and it was even shown that a continuous time quan- 
tum walk on a specific graph can be used for exponential 
algorithmic speed-up 

All studies on quantum walks so far have, however, 
been based on a single walker. In this article we study a 
discrete time quantum walk on a line with two particles. 
Classically, random walks with K particles are equivalent 
to K independent single-particle random walks. In the 
quantum case though, a walk with K particles may con- 
tain entanglement, thus offering a resource unavailable 
in the classical scenario which can present interesting ad- 
vantages. Moreover, in the case of identical particles we 
have to take into account the effects of quantum statis- 
tics, giving an additional feature to quantum walks that 
can also be exploited. In this work we explicitly show 
that a quantum walk with two particles can indeed be 
tuned to behave very differently from two independent 
single-particle quantum walks. This paves the way for 
new quantum algorithms based on richer quantum walks. 



Let us start by introducing the discrete time quantum 
walk on a line for a single particle. The relevant degrees 
of freedom are the particle's position i (with igZ) on 
the line, as well as its coin state. The total Hilbert space 
is given by 7i = Tip <g> Tic, where Tip is spanned by 
the orthonormal vectors representing the position 

of the particle, and Tic is the two-dimensional coin space 
spanned by two orthonormal vectors which we denote as 
IT) and ||). 

Each step of the quantum walk is given by two sub- 
sequent operations. First, the coin operation, given by 
Uc S SU(2) and acting only on Tic, which will be the 
quantum equivalent of randomly choosing which way the 
particle will move (like tossing a coin in the classical 
case) . The non-classical character of the quantum walk 
is precisely here, as this operation allows for superposi- 
tions of different alternatives, leading to different moves. 
Then, the shift-position operation S will move the par- 
ticle accordingly, transferring this way the quantum su- 
perposition to the total state in Tt. The evolution of the 
system at each step of the walk can then be described by 
the total unitary operator: 



u = s(i P ® u c ), 



(i) 



where Ip is the identity operator on Tip. Note that if a 
measurement is performed after each step, we will revert 
to the classical random walk. 

In this article we choose to study a quantum walk with 
a Hadamard coin, i.e. where Uc is the Hadarmard oper- 
ator H: 



1 1 
1 -1 



(2) 
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Note that this represents a balanced coin, i.e. there is a 
fifty-fifty chance for each alternative. The shift-position 
operator is given by: 

S = (£ |» + 1)<^ ® ITXTI + fe I* - Wl) ® 

^ ' ' . (3) 

Therefore, if the initial state of our particle is, for in- 
stance |0) ® |t), the first step of the quantum walk will 
be as follows: 

|o>®|t) |o)®i(|t) + ||)) 

A ^(|l)®|t) + |-l)®|i)). (4) 

We see that there is a probability of 1/2 to find the par- 
ticle in position 1, as well as to find it in position 2, just 
like in the classical case. Yet, if we let this quantum walk 
evolve beyond the (two) initial steps before we perform a 
position measurement, we will find a very different prob- 
ability distribution for the position of the particle when 
compared to the classical random walk, as it can be seen 
in Fig. [2 for N = 100 steps. 



earlier for the one-particle case. Note also that in the 
case of identical particles we have to restrict H12 to its 
symmetrical and antisymmetrical subspaces, respectively 
for bosons and fermions. 

Let us then consider the case where both particles start 
the quantum walk in the same position, 0, but with differ- 
ent coin states ||) and \\). In the case where the particles 
are in a pure separable state, our system's initial state 
will be given by: 



Mf)i2 = |0,I)i|0,T); 



(6) 



We can also consider an initial pure state entangled in the 
coin degrees of freedom. In particular, we will consider 
the following two maximally entangled states: 

= -L(|0,i>i|0,t)2±|0,t)i|0,i)a), (7) 

differing only by a relative phase. Note that if we were 
considering identical particles on the same point, our 
system would have to be described by these states, for 
bosons and fermions respectively. 

Each step of this two-particles quantum walk will be 
given by: 
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FIG. 1: Probability distribution for a classical random walk 
(a) on a line after N = 100 steps, as well as for a quantum 
walk (b) with initial state |0) g) |T) and a Hadamard coin. 

Let us consider the previous quantum walk, but now 
with two non-interacting particles on a line (not necessar- 
ily the same). If the particles are distinguishable and in 
a separable state, the position measurement of one par- 
ticle will not change the probability distribution of the 
other, they are completely uncorrelated. On the other 
hand, if the particles are entangled, a new resource with 
no classical equivalent will be at our disposal. 

The joint Hilbert space of our composite system is 
given by: 

H12 = Hi®H 2 = {Hpa <g> H c ,i) ® (Wp,2 ® Hc,2), (5) 

where 7i\ and TL2 represent the Hilbert spaces of particles 
1 and 2 respectively. Since the relevant degrees of free- 
dom in our problem are the same for both particles, we 
have that both Tii and TI2 are isomorphic to TL defined 



(8) 



where U is given by equation and is the same for both 
particles. After N steps, the state of the system will be, 
in the case of the initial conditions ©: 



I^)i2 



^1^)12 = ^10,1)1^10,1)2. (9) 



Fig. |2]a shows the joint probability distribution 
Pf 2 (i,j; N) for finding particle 1 in position i and particle 
2 in position j for N — 30 steps. Note that, since the par- 
ticles are uncorrelated, Pf 2 (hj', N) is simply the product 
of the two independent one-particle distributions: 



Pg(i,j;N) = P?(i;N) x P 2 s (j;N), 



(10) 



where Pf (i; N) is the probability distribution for finding 
particle 1 in position i after N steps, and similarly for 
P2 (j;N) and particle 2. This can also be observed in 
Fig.|5Ja, which is clearly the product of two distributions 
like the one in Fig. ^b, one biased to the left for particle 
1 and the other to the right for particle 2, accordingly 
with the initial conditions given by equation ©. 

In the case of entangled particles, the state of the sys- 
tem after N steps will be: 



'12 



(11) 



I^)i2 - u 12 \^ 
-L (u"\Q, [) 1 T ) 2 ± LT^IO, T)i U N \0, |) 2 



The probability distribution for finding particle 1 in po- 
sition i and particle 2 in position j in the "+" case, 
P^ 2 (i,j;N), is represented in Fig.EJb, for N = 30. Sim- 
ilarly, Fig. [21c shows the distribution P^(i,j;N) of the 
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"-" case, again for N = 30. The effects of the entangle- 
ment are striking when comparing the three distributions 
in Fig. In particular, we see that these effects signif- 
icantly increase the probability of the particles reaching 
certain configurations on the line, which otherwise would 
be very unlikely to be occupied. In all cases, the max- 
ima of the distributions occur around positions ±20. In 
the "+" case it is most likely to find both particles to- 
gether, whereas for "-" the former situation is impossible, 
and the particles will tend to finish as distant as possible 
from one another. 

Let us now consider the individual particles in the 
entangled system. The state of each of them can be 
described by the reduced density operator pi, 2 (N) = 
Tr 2 .\ (|'0w)i2 12(V'wI)j which consists of an equal mix- 
ture of the one-particle states U N \0,l) and U N \0, f)- 
Thus, given one particle, the probability to find it in po- 
sition i after N steps is given by the following marginal 
probability distribution: 

P 1 ± S;N) = ±[P l (z;N) + P ] (t;N)]. (12) 

where P±(i;N) is the probability distribution for find- 
ing the particle in state U N \Q, J.) in position i after TV 
steps, and similarly for P^(i;N) and a particle in state 
U N \0,1). Note that in the case of initial conditions © 
we have the marginal probabilities Pf(i;N) = P±(i;N) 
and P 2 (i; N) = P-\{i; N). But now, contrary to the sep- 
arable case, the joint probability P^(i,j; N) is no longer 
the simple product of the two one-particle probabilities, 
as it contains information about the non-trivial corre- 
lations between the outcomes of the position measure- 
ment performed on each particle. Therefore, to investi- 
gate more quantitatively the difference between quantum 
walks with two distinguishable particles in a pure separa- 
ble state and two entangled particles in their coin degree 
of freedom, we must look at joint (two-particle) rather 
than individual properties. 

First, let us consider the distance between the two par- 

S i 

tides, . For that, we define the random variables Xi 
and x 2 as the outcomes of the (single particle) measure- 
ment of the position of particle 1 and particle 2 respec- 
tively, after N steps of the quantum walk (they can take 
integer values between N and —N). We can now define 
the distance as, for the three different initial conditions: 

Af^ = \xi-x 2 \. (13) 

The expectation value (Af 2 ) of the distance is given in 
Table for different N. We see that in the "-" case the 
particles tend, on average, to end the quantum walk more 
distant to each other, whereas in the "+" case they tend 
to stay closer, and somewhere in between in the separable 
case. In fact, for fixed N, we always have: 

(Ar 2 )-<Af 2 ) = (Af 2 )-<A+). (14) 

Let us now consider the correlation function between 




A. Separable initial conditions 




C. Entangled initial conditions "-" 



FIG. 2: Two-particles probability distributions after N = 
30 steps for different initial conditions: (A) separable state 
|i/>o)i2; (B) |V'o")i2 state; and (C) |$f)i2 state. Note the 
different vertical ranges. 

the spatial distribution of each of the two particles: 

C s - ± (x 1 ,x 2 ) = (x lX2 ) - (x 1 )(x 2 ). (15) 

Clearly, in the case of the separable state © this corre- 
lation is always zero. For maximally entangled particles, 
the values of C ± {x\,x 2 ) are presented in Table HT1 for dif- 
ferent N. Given the symmetry of the P-f^ distributions, 
we have that in those cases (x\) — (x 2 ) = 0. Thus, the 
sign difference in the correlation function expresses the 
tendency for the two particles in the "-" case to end the 
quantum walk on different sides of the line (with respect 
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Expectation value (Af^) after TV steps 



Nb 


of steps N 


10 


20 


30 


40 


60 


100 


Init. 


cond. \iPq)\2 


8.8 


17.5 


26.0 


34.9 


52.2 


87.0 


Init. 


cond. |V>o)i2 


7.1 


14.7 


21.9 


29.5 


44.3 


73.9 


Init. 


cond. \i/)q)i2 


5.5 


11.9 


17.8 


24.1 


36.3 


60.8 



TABLE I: Average distance (Afp) after N steps. 



Correlation function C (xi, X2) after N steps 



Nb. 


of steps N 


10 


20 


30 


40 


60 


100 


Init 




-16.8 


-69.8 


-153.5 


-276.2 


-619.7 


-1718.3 


Init 


C |V>0 S )l2 




















Init 


C- |^o")l2 


4.8 


7.3 


13.7 


15.1 


23.1 


39.1 



TABLE II: Correlation function C s,± (a;i, X2) after N steps. 



to its origin 0), and on the same side for the "+" case. 

Finally, let us now calculate, for the different initial 
conditions, the probability of finding at least one parti- 
cle in position i after N steps: 'P S ' ± (i; N). This is clearly 
a joint property as it depends on both one-particle out- 
comes: 



V b ' ± (i;N) 

N 



7=-JV 



= [if' ± (i;JV) + if V) 
= [PMW+P&N)] 



(16) 



Given a one-particle probability distribution, say 
Pl(i;N) [note that P^(i;N) = P L (-i;N)], our prob- 
ability V s,± {i] N) decreases with the joint probability 
P 12 ' (i,i;N) and is maximal in the "-" case, as we al- 
ways have P^{i,i;N) — 0. In fact, around the points 
(20, 20) and (-20, -20) in Fig. Hwe clearly have: 



V-(i;N) > V s (i;N) > V + (i;N). 



(17) 



We see that, by introducing entanglement in the initial 
conditions of our two-particles quantum walk, the 
probability of finding at least one particle in a particular 
position on the line can actually be better or worst than 
in the case where the two particles are independent. 
Note that in this case this does not depend on the 
particular amount of entanglement introduced, as both 
states in Eq. Q are maximally entangled, but rather on 
their symmetry/relative phase. 

This work opens way for generalization in a number of 
ways. First, one could consider periodic or other bound- 
ary conditions on the line, or more general graphs 
Note that the positions on the line of our two particles 
could also be interpreted as the position of a single parti- 
cle doing the quantum walk on a regular two-dimensional 



lattice. More general coins could also be considered 0,0] , 
including entangling and non-balanced coins, as well as 
different initial states. One could also augment the num- 
ber of particles, study these quantum walks in continu- 
ous time or in their asymptotic limit. Furthermore, also 
very interesting and promising is to investigate the use 
of multiparticle quantum walks in the design of quantum 
algorithms 8], or in solving mathematical or practical 
problems that could be encoded has a quantum walk, 
such as the estimation of the volume of a convex body 
or the connectivity in a P2P network [To| . 

Finally, some brief comments about implementations 
of our two-particles quantum walk on a line. The meth- 
ods recently proposed for the sin gle-p article case usin 
cavity QED [ll|, optical lattices [12( or ion traps [i 
could be adapted to our two-particles case. For instance, 
in the latter we could encode the coin states in the elec- 
tronic levels of two ions and the position in their com or 
stretch motional modes: the coin flipping could then be 
obtained with a ir/2 Raman pulse and the shift with a 
conditional optical dipole force Another possibility 
is to send two photons through a tree of balanced beam 
splitters which implement both the coin flipping and the 
conditional shift, aga in g eneralizing a scheme proposed 
for a single particle |6j,ll5|. Note that this could be imple- 
mented with other particles as well, e.g. electrons, using a 
device equivalent to a beam splitter [16| . Also very inter- 
esting, regardless of any particular method or technology, 
is the possibility of using two indistinguishable particles 
on the same line to implement our quantum walk. Say 
we encode the coin degrees of freedom in the polarization 
of two photons, or in the spin of two electrons: if the two 
particles start in the same position 0, then they will be 
forced to be in state given by Eq. J7J ( "+" for bosons and 
"-" for fermions). Although the particles will initially be 
only entangled in the mathematical sense (as they cannot 
be addressed to extract quantum correlations), this is a 
perfectly valid way of preparing the |'0o")i2 and |V'r7)i2 
initial states for our quantum walk, saving us the trouble 
of generating entangled pairs. Thus, the indistinguisha- 
bility of identical particles appears once again as a re- 
source for quantum information processing [TtI ]. here in 
particular offering a way to simplify the preparation of 
the initial states for a two-particle quantum walk on a 
line. 

In this article we introduced the concept of a quantum 
walk with two particles and studied it for the case of a dis- 
crete time walk on a line. Having more than one particle, 
we could now add a new feature to the walk: entangle- 
ment between the particles. In particular, we considered 
initial states maximally entangled in the coin degrees of 
freedom and with opposite symmetries, and compared 
them to the case where the two particles were initially 
unentangled and thus independent. We found that the 
entanglement in the coin states introduced spatial cor- 
relations between the particles, and that their average 
distance is larger in the "-" case than in the separable 
case, and is smaller in the "+" case. This could benefit 
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algorithmic applications which require two marked sites 
to be reached which are known a priori to be on the oppo- 
site or on the same sides of a line. We also found that the 
introduction of entanglement could increase or decrease 
the probability to find at least one particle on a given 
point of the line. This increase could allow us to reach 
a marked site faster than with two unentangled quan- 
tum walkers. The entanglement in the initial conditions 
thus appears as a resource that we can tune according to 
our needs to enhance a given application (algorithmic or 
other) based on a quantum walk. 
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